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ABSTRACT 

The combination of several mass-decreasing processes may critically affect the struc- 
ture of open clusters (OCs), to the point that most dissolve into the field in a time-scale 
shorter than w 1 Gyr. Therefore, as is observed, old OCs should be sparse within the 
Galaxy. These mass-decreasing processes can only be quantifiably examined given 
the derivation of the fundamental parameters, age, distance and mass, for a sample 
of old OCs. In this paper we focus on 4 candidate old star clusters, namely, LKl, 
LKIO, FSR1521, and FSR1555. The first two clusters, discovered by Le Duigou & 
Knodlseder, are projected towards the Cygnus Association, whilst the remaining two 
have been detected in the 4th quadrant as cluster candidates from stellar overdensi- 
ties by Froebrich, Scholz & Raftery. To analyse the target clusters we construct near 
infrared colour-magnitude diagrams (CMDs) and derive stellar radial density profiles 
(RDPs). The CMDs are constructed using 2MASS J, H, and Ks bands, and the in- 
trinsic morphologies of the target OCs within these diagrams are revealed by applying 
a field-star decontamination algorithm. Fundamental parameters are estimated with 
Padova isochrones built for the 2MASS filters. We derive extinctions to the objects 
within the range 3.4 ^ Ay ^ 8.9, which makes them suitable for the near-infrared 
analysis, ages within 1.0 — 2.0 Gyr, and distances from the Sun within 1.4 — 4.5 kpc. 
These distances, in conjunction with the positions in the sky, place the present 4 
OCs close to the solar circle (< 0.6 kpc). For LKIO our photometry reaches a depth 
« 3 mags below the main sequence turn off, from which we derive a relatively steep 
mass function slope (x — 2.4 ± 0.4) when compared to the Salpeter value (x = 1.35). 
LK 10 is a rather massive old OC, with a mass within 1360 ^ m{MQ) ^ 4400, for 
stars in the observed magnitude range and the extrapolation to 0.08 Mq, respectively. 
The mass estimated in the restricted magnitude range for the remaining, more distant 
OCs is within 260 ^ m{MQ) ^ 380. However, similarity with the CMD morphology 
and red clump of LK 10 suggests that they may be as massive as LK 10. The RDPs 
are well represented by a King-like function, except LK 10, which has a central cusp 
suggesting post-core collapse. Structurally, LKl, FSR1521, and FSR155 are similar 
to a sample of nearby OCs of comparable age. 

Key words: ( Galaxy:) open clusters and associations: general; ( Galaxy:) open clus- 
ters and associations: individuahLK 1, LKIO, FSR1521, and FSR1555. 



1 INTRODUCTION 

The majority of the Galactic open clusters (OCs), the less 
massive ones in pa r ticular, do not survive t he k, 1 Gyr age 
barrier (|Friellll995l : iBonatto fc Bical l2007bl . and references 
therein). These systems are affected by several processes 
that continually change the gravitational potential, such as 



* E-mail: charles@if.ufrgs.br 
f E-mail: bica@if.ufrgs.br 



mass loss associated with stellar evolution, mass segrega- 
tion and evaporation, tidal interactions with the Galactic 
disc and bulge, and collisions with giant molecular clouds. 
As clusters age, these mechanisms accelerate the internal 
dynamical evolution, which leads to important changes in 
the structure and flat, eroded mass functions. In most 
cases, OCs have their stellar content completely dissolved in 
the Galactic field, or leave only poorly-populated remnants 
(|Pavani fc Bicair2007l and references therein). 

Ample evidence gathered on theoretical (e.g. 
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Spitzerl [l95i: iLamers fc Gieled |2006|) 



„ ------ | .- N-body (e.g . 

Baumgardt fc Makind l2003l: iGoodwin fc BastianI l20od : 
Khalisi. Amar o- S coanc fc Spurzcm' '2007'). and observa- 
tiona l (e.g. [v an den B crgh 1957; Oort 19581 : Ivon Hoernej 
1 19581 : iPiskunov et"all |2007| ) grounds indicate that, near 
the Solar circle, the disruption-time scale (tdis) is shorter 
than ~ 1 Gyr. The disruption-time scale around the 
Solar circle a ppears to depend on mass as idis ~ M°'^^ 
ijLamers fc Gieles,2006 ). Thus, for clusters with mass within 
10^ - IO^Mq, a disruption time of 75 < taisiMyr) < 300 
should be expected. In addition, disruption processes are 
stro nger for the more centr ally located and lower-mass OCs 
fsee lBonatto fc BicallioOTbl for a review on disruption effects 
and time-scales). Gyr-class OCs are found pr eferentially 
near the Solar circle a nd in the outer Galaxy (e.g. lFrielll995l : 
iBonatto et al.l [2006al ). where the frequency of potentially 
disrupting dynamical interactio ns with giant molecular 
clouds and the disc is lower (e.g. ISalaris. Weiss fc Percivall 
l2004l : lUpgrenT Mesrobian fc Kerridgdligm " As an extreme 
case, inner (-Rgc ^ 150 pc) Galactic tidal fields can dissolve 
a ma ssive st ar cluster in a t ime-scale as short as ~ 50 Myr 
ifPort egi es Zwart et al.ll2002l ) . 

What should be expected from the above scenario 
is that only a small fraction of the OCs survive the 
Gyr age-barrier, with the successful ones spending most 
of their existences preferentially at large Galactocentric 
distances. Indeed, present-day statistics show that of the 
~ 1000 OCs with known age listed in the webdaQ 
database, 180 are older than 1 Gy r, and only 18 (~ 
2%) are older than 4 Gyr (se e also lOrtolani et al.1 l2005l : 
lOrtolani. Bica fc Barbuvl l2005h . Additionally, most of the 
OCs older than 1 Gyr so far identified are located outside 
the Solar circle (see, e.g., the spatial distribution of OCs of 
different ages in Fig. 1 of lBonatto fc Bicall2007bh . 

Given the aspects discussed above, a more thorough 
and wide-ranging exploration of the OC parameter space - 
the older ones in particular - is fundamental for a better 
understanding of the dynamical process that affect the star 
clusters. In this sense, it is important to unambiguously es- 
tablish the nature of as many old OCs as possible, followed 
by the derivation of reliable astrophysical parameters. Such 
parameters, in turn, can be used in studies of star formation 
and evolution processes, dynamics and survival rates of star 
clusters, the geometry of the disc, among others. 

In the present paper we focus on 4 candidate old OCs, 
namely, LKl, LK 10, FSR1521, and FSR1555. LK 1 and 
LK 10 are projected towards C ygnus and were found by 
iLe Duigou fc Knodlsedeij (|2002h . An additional interesting 
point is to determine whether LK 1 and/or LK 10 are part of 
the Cygnus Association. FSR 1521 and FSR 1555 were found 
as candidate open clusters inferred fro m stellar overdensities 
bv lFroebrich. Scholz fc Raftervl (|2007l ). 

This work employs near-IR J, H, and Ks photometry 
obtained from the 2MAS^| Point Source Catalogue (PSC). 
The spatial and photometric uniformity of 2MASS, which 
allow extraction of large surrounding fields that provide 
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Table 1. Previous identifications of old clusters by our group 



Cluster 


Age 

(^yr) 


Reference 




FSR 869 


1.5 


Bonatto & Bica (20080) 




FSR 942 


1.0 


Bonatto & Bica (20080) 




FSR 70 


>5 


Bica. Bonatto &; Camareo 


[2008) 


FSR 124 


1.0 


Bica. Bonatto &; Camareo 


[2008) 


ES0 277SC1 (FSR 1723) 


0.8 


Bica, Bonatto & Camarso (2008) 


FSR 1737 


>5 


Bica, Bonatto & Camargo 


f2008) 


RulOl (FSR 1603) 


1.0 


Bica & Bonatto ('2008) 




FSR 1744 


1.0 


Bonatto & Bica ('2007b) 




FSR 89 


1.0 


Bonatto & Bica ('2007b) 




FSR 31 


1.1 


Bonatto & Bica (■2007b) 




Cz23 (FSR 834) 


4.5 


Bonatto & Bica f2008d) 




FSR 1716 


>5 


Bonatto & Bica f2008d) 




FSR 1415 


2.0 


Momanv et al. f2008) 




Bica 6 (BBSl) 


1.0 


Bonatto. Bica &; Santos Jr. 


C2008) 








Table Notes. Bica6 (DAML02), located at «(2000) = 
05'"06™20* and <5(2000) = +39°09'50", is probably physi- 
cally associated with the planetary nebula PK 167-0.1. 



high star-count statistics, make it an excellent resource to 
gather photometric data on a broad variety of star clus- 
ters, the wide field ones in particular. For this purpose we 
have developed quantitative tools to statistically disentan- 
gle cluster evolutionary sequences from field stars in colour- 
magnitude diagrams (CMDs), which are subsequently used 
to investigate the nature of star cluster candidates and to 
deriv e astrophysical parameters of t he confirmed clusters 
(e.g. iBica. Bonatto fc Camargol boOSi ) . Basically, we apply 
(i) field-star decontamination to quantify the statistical sig- 
nificance of the CMD morphology, which is fundamental to 
deriving reddening, age, and distance from the Sun, and 
(a) colour-magnitude filters, which are essential for intrin- 
sic stellar radial density profiles (RDPs), as well as lumi- 
nosity and mass functions (MPs). In particular, the use of 
field-star decontamination in the construction of CMDs has 
proved to constrain the age and distance more than when 
working with the raw (observed) photo metry, especially for 
low-latitude OCs (|Bonatto et al.ll2006al ). 

This paper is organised as follows. In Sect. [2] we recall 
recent additions to the known old OCs made by our group. 
Sect. |3]contains basic properties and reviews literature data 
(where available) on the present star cluster candidates. In 
Sect.Uwe present the 2MASS photometry, build CMDs, and 
apply the field-star decontamination algorithm. In Sect. [5] 
we derive cluster fundamental parameters. Sect. [6]describes 
cluster structure by means of stellar RDPs. In Sect. [7] we 
provide estimates of cluster mass. In Sect.|8]we compare the 
structural parameters and dynamical state of the present 
clusters with those of a sample of nearby OCs, we also dis- 
cuss effects of the location in the Galaxy on their structure. 
Concluding remarks are given in Sect. (9) 
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Figure 1. Top panels: 5' X 5' 2MASS Ks images of LK 1 (left) and LK 10 (right). Bottom: same for FSR 1521 (loft) and FSR1555 
(right). Images provided by the 2MASS Image Service. The small circle indicates the central coordinates (cols. 4 and 5 of Table [2}. 
Figure orientation: North to the top and East to the left. 



2 RECENT ADDITIONS TO THE KNOWN 
POPULATION OF OLD OCS 

In recent years our group has been systematically analysing 
infrared clusters or candidates, establishing their nature and 
deriving cluste r fundamental parameters with the 2MASS 
catalogue (e.g. iBica. Bonatto fc Camargol f2008l ) . We make 



use of a field decontamination algorithm (described in 
Sect. 14.111 to statistically extract estimated cluster sequences 
from CMDs. Especially for crowded fields, the cluster se- 
que nce isolation requires some form of membership selection 
(e.g. lBonatto fc Bica,2007bi ). Old OCs are intrinsically rarer 
than young clusters (Sect, [ij, but the 2MASS catalogue, 
coupled to the statistical tools that we have developed to 
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Table 2. Previous data and present results on the clusters 



Literature This paper 



Cluster 


a(2000} 


5(2000) 


q(2000) 


S(2000) 


f 


h 


Age 


Av 


dQ 


«GC 




''GC 


^GC 




(hms) 


(°"') 


(hma) 


(° "') 


(°) 




(Gyr) 


(mag) 


(kpe) 


(kpc) 


(kpe) 


(kpc) 


(kpc) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(S) 


(£>) 


(10) 


(11) 


(12) 


(13) 


(14) 



LK 1 20:24:25 +36:24:18 

LKIO 20:39:13 +39:58:37 

FSR1521 09:55:23 -56:36:06 

FSH1555 10:49:07 -59:03:17 



(t) +36:24:30.0 75. 

(t) +39:58:43.2 79. 

09:55:22.3 -56:37:16.8 280 

10:48:55.8 -59:02:58.9 287 



24 -0.69 1.0 ± 0.2 8.9 ± 

84 -0.92 1.0 + 0.1 8.5 + 

.45 -1.64 2.0 + 0.5 3.4 + 

.75 +0.16 1.5 + 0.5 3.7 + 



3 4.0 + 0.4 7.3 + 0.2 -6.2 

3 1.4 + 0.1 7.1 + 0.1 -7.0 

5 4.5 + 0.7 7.8 + 0.4 -6.4 

3 4.1 + 0.2 7.1 + 0.1 -6.0 



+ 0.1 3.9 + 0.4 -0.05 + 0.01 

+ 0.1 1.4 + 0.1 -0.02 + 0.01 

+ 0.1 -4.4 + 0.7 -0.13 + 0.02 

+ 0.1 -3.9 + 0.2 0.01 + 0.01 



Table Notes. Coordinates (eols. 2 and 3) of LK 1 and LKIO from fLe~Duigou Knodlse^el 120021 ): FSR1521 and FSR 1555 from iFroebrich. Scholz Raftervl 120071) : 
Cols. 4-7: optimised coordinates; (f): same value as in th e literature; Col. 9: reddening towards the cluster's central region (Sect.lsT. Col. 10: distance from the Sun. 
Col. 11: cluster Galactocentric distance for Kq = 7.2 kpc IBica et £101200(3 ). Cols. 12-14: coordinate components projected onto the Galactic plane. 



deal with clusters and comparison fields allow the detection 
(and deri vation of reliable astrophysical parameters) of o lder 
OCs (e.g. iBonatto fc Bi ca"2007b'; 'Bica fc Bonattol lBoj ). 

Recently, IFroebrich . Scholz & Raftory (20q3) provided 
a catalogue of star cluster candidates corresponding to 1021 
stellar overdensities detected in the 2MASS database. This 
catalogue covers \b\ < 20° and all Galactic longitudes, 
and has become an important source of new star clus- 
ters. Several foUowup studies have explored the FSR cat- 
alogue with different appr oaches, revealing new globular 
clust ers, such as FSR 173 5 iFroebrich. Meusinger fc Schola 
l2007l ) and FSR 1767 (|Bonatto et al.l l2007f) . and the 
probable GCs FSR 584 jB ica ct al' |2007| ) and FSR 190 
IFroebrich. Meusinger fc Davis 2007. ) . We show in Table [1] 
14 recent identifications of old OCs made by our group. 
Three of them are optical objects, while the remaining 11 
are infrared ones. Compared to the DAMLO^ optical cata- 
logue of OCs, our results increased the known sample of old 
OCs by « 5%. 

In the present paper we merge the old OCs in Table [1] 
with the classical optical old OCs (WEBDA) for comparison 
purposes with the 4 objects dealt with in this paper. This 
analysis is given in Sect. [8] 



(|2007l ) as a highly probable star cluster candidate. They 
derived the core and tidal radii (measured in 2MASS H 
images) = 1^6 and J?f = 25!3, respectively. FSR 1521 
can be seen in the 5' x 5' Ks image shown in Fig. [1] (bottom- 
left panel). 

FSR 1555 was als o classified as a highly proba b le sta r 
cluster candidate by IFroebrich. Scholz fc Raftervl l|2007h , 
who derived = 1^7 and R^ = 11.'6 for this object. 
FSR 1555 is shown in the 5' x 5' Ks image shown in Fig. [T] 
(bottom-right panel). In general, the clusters in the present 
sample are very contaminated by field stars, which requires 
specific tools to analyse them. 

Table [2] provides fundamental data on the objects, 
where the literature coordinates are given in cols. 2 and 
3. However, when we built the RDPs based on these co- 
ordinates (Sect. we noticed that, in all cases, the coor- 
dinates where the maximum stellar number-density occurs 
are slightly shifted with respect to the literature positions. 
Thus, hereafter we will refer as cluster coordinates those that 
maximise the central stellar density (given in cols. 4-7). The 
age, central reddening, distance from the Sun, Galactocen- 
tric distance, and the components projected onto the Galac- 
tic plane derived in the present study (Sect. [5]) are given in 
Cols. 8 to 14. 



3 THE TARGET CLUSTERS 

iLe Duigou fc Knodlsedeil l|2002l ) provided a list of 17 star 
clusters and candidates in the Cygnus direction. As they 
point out, 12 of these object s (hereafter de s ignat ed by LK) 
had already been found by iDutra fc Bical (1200 ll ). Most of 
the LK objects appear to be embedded clusters, but some 
of them, like LK 1 and LK 10, have features typical of older 
clusters. With 2MASS photometry, they were able to es- 
timate some cluster parameters. For LK 1 they derived a 
radius (containing 90% of the stars) Rgo = 3.'0, the dis- 
tance modulus DM — 11.0, the absorption in the K band 
Ak = 1.0 - 2.0, and a mass within M = 1100 - 3130 M©; 
they suggested that LK 1 may be a rather evolved clus- 
ter. As for LKIO, they found R90 = 6^7, DM = 11.0, 
Ak = 0.6 - 1.2, M = 1010 - 4100 M©, a rather steep mass 
function (MF: <;/>(m) tx m"'^+^') slope x = 1-98 ± 0.25, and 
suggested that it is quite evolved. Near-infrared 2MASS Ks 
images of LK 1 and LK 10, covering 5' x 5' fields, are shown 

in Fig. [1] (top panels) . 

FSR 1521 was classified by IFroebrich. Scholz fc Raftervl 



4 COLOUR-MAGNITUDE DIAGRAMS WITH 
2MASS PHOTOMETRY 

Photometry in the 2MASS J, H, and Ks bands was ex- 
tracted in circular fields of radius J?cxt centred on the co- 
ordinates of the objects (Table [2)| by means of VizieE[f|. 
J?oxt should be large enough to allow the determination 
of the background level (Sect. [S}. In the present cases, 
R^^^ = 30' (FSR1521 and FSR1555), .Rcxt = 40' (LKl), 
and iicxt = 60' (LKIO), which are considerably larger than 
the cluster radius (Sect. |S] and col. 5 of Table |3|. Previ- 
ous works with OCs in different environments (Sect.[l| have 
shown that in the absence of a populous neighbouring clus- 
ter and significant differential absorption (Sect. H7l.l|l . wide 
extraction areas provide the necessary statistics for a con- 
sistent colour and magnitude characterisation of the field 
stars. For decontamination purposes, comparison fields were 
extracted within wide rings located beyond the cluster radii. 
As photometric quality constraint, the 2MASS extractions 
were restricted to stars (i) brighter than those of the 99.9% 



^ http://www.astro.iag.usp.br/ wilton/ 



I http: //vizier, u-strasbg.fr/ viz- bin/VizieR ?-source= II/246\ 
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LK1:Raw LK10:Raw 




Figure 2. Stellar surface-density (T{stars arcmin"^) of LK 1 (left 
panels) and LK 10 (right). The curves were computed for a mesh 
size of 3' X 3', centred on the coordinates in Table[2l The observed 
(raw) and field-star decontaminated photometry are shown in the 
top and bottom panels, respectively. 

Point Source Catalogue completeness limil[f| in the cluster 
direction, and (ii) with errors in J, H, and Ks smaller than 
0.3 mag. The 99.9% completeness limits refer to field stars, 
and depend on Galactic coordinates. For the present clus- 
ters, the fraction of stars with J, H, and Ks uncertainties 
smaller than 0.06 mag is « 80%. A typical distribution of 
uncertainties as a function of magnitude, for clusters pro- 
je cted towards the c entral parts of the Galaxy, can be found 
in lBonatto fc Bical (|2007al ). Reddening transformations use 
the relations Aj/Ay = 0.276, Ah I Ay = 0.176, Ak^/Av = 
0.118 , and Aj = 2.76 x E{J - H) (|Dutra. Santiago fc Bical 
|2002| ). for a constant total-to-selective absorption ratio 
Rv = 3.1. These ratios were derived from the extinction 
curve of lCardeUi. Clayton fc Mathid l|l989D . 

CMDs displaying the J x [J - H) and J x {J ~ Ks) 
colours built with the raw photometry of the present clus- 
ters are shown in Figs.|3]-[7](top panels). For all clusters, the 
sampled region is larger than the respective core (Table U). 
When qualitatively compared with the CMDs extracted 
from the equal-area comparison fields (middle panels), fea- 
tures typical of old OCs are apparent. A relatively populous 
red clump (at 13.8 < J < 14.7 and 1.3 < (J ~ H) < 1.5) and 
a main-sequence turn off (MSTO) stand out over the field 
contamination of LK 1 (Fig. A somewhat less-populous 
red clump (11.3 < J < 12.7, 1.2 < {J - H) < 1.5) and 
about 3 MS mags (J > 13.4) are seen in LK 10 (Fig.[Sl). LK 1 
(Fig. 11), FSR1521 (Fig.EJ, and FSR1555 (Fig. present 
similar CMDs, with clear red clumps, the MSTO, and about 
1 mag of the MS below. The red clump of FSR 1521 occurs 



^ According to the 2MASS Level 1 Requirement, at 
://www.ipac.caltech.edu/2mass/releases/allsky/doc/\ 



FSR1555:Raw FSR1521:Raw 




Figure 3. Same as Fig. [2] for FSR 1555 (left) and FSR 1521 
(right). 



at 12.3 < J < 13.2 and 0.8 < {J - H) < 1.0, while for 
FSR 1555 it is at 12.3 < J < 13.5 and 0.7 < {J - H) < 1.0. 

Finally, we show in Figs. [5] and |3] the spatial distri- 
bution of the stellar surface-density as measured in the 
2MASS photometry. We compute the surface density a (in 
units of stars arcmin"^), in a rectangular mesh with cells 
of dimensions 3' x 3'. The meshes reach total offsets of 
Aa[ = \AS\ « 20' with respect to the centre (Table 
in right ascension and declination; for LK 10 we use offsets 
of 30'. In all clusters, the core (Table |4} is contained in the 
central cell. 

With respect to the surface-densities built with the ob- 
served (raw) photometry (top panels of Figs. [2] and |3} , an 
important excess appears in the central cell, except for LK 10 
which, because of the contamination by disc stars, presents a 
rather irregular distribution. FSR 1521 and FSR 1555, on the 
other hand, clearly detach in the central cell (Fig. [Sjl against 
more uniform surrounding fields. As shown in the bottom 
panels, the cluster overdensities are clearly enhanced with 
respect to the surroundings in the field-star decontaminated 
surfaces (Sect. HTTjl . 

4.1 Field-star decontamination 

As expected of low- latitude clusters (Table the stellar 
surface-density in the direction of the objects (Figs. [2] and [3]) 
clearly shows that field-star contamination, essentially from 
disc stars, should be taken into account. This fact is con- 
firmed by the qualitative comparison between the CMDs ex- 
tracted within the cluster and in the field (Figs. |4][7]). Thus, 
the field-star contribution should be quantified in each case 
to better define the intrinsic CMD morphology. 

Field-star decontamination is a very important, yet dif- 
ficult, step in the identification and characterisation of star 
clust ers. Several approa ches have been used to this purpose 
(e.g. lMercer et allbOOSf ). and most of them are based essen- 
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Figure 4. 2MASS CMDs extracted from the R < 2' region 
of LK 1. Top panels: observed pliotometry with the colours 
J X {,] ~ H) (left) and J X (J — .ft's) (right). Middle: equal-area 
(29.93' < R < 30') extraction from the comparison field, where 
the disc contamination is present. Bottom panels: decontaminated 
CMDs that suggest a relatively reddened and distant MSTO, red 
clump, and giant branch typical of old OCs, fitted with the 1 Gyr 
Solar-metallicity Padova isochrone. The shaded polygon corre- 
sponds to the colour-magnitude filter (Sect. |6ll. Arrows in the 
bottom panels show the reddening vector computed for Ay = 2. 



tially on two different premises. The first works with spatial 
variations of the star-count density, but does not take into 
account CMD properties. In the latter, stars in a CMD ex- 
tracted from an assumed cluster region are subtracted ac- 
cording to colour and magnitude similarity with the stars 
of an equal-area comparison field CMD. These methods, to- 
gether with the one we work with, are based on photomet- 
ric properties only. Ideally, more robust results on cluster 
membership determination would be obtained if another in- 
dependent parameter, such as the proper motion of member 
and comparison field stars, is taken into account. However, 
for proper motions t o be useful the target cluster s hould be 
relatively close (e.g. lAlessi. Moitinho fc Diasllioosi ) and/or 
to have been observed in widely-apart epochs, preferentially 
with high resoluti on, such as in the case of the globular 
cluster NGC6397 (jRicher et al.ll2008h . Neither condition is 
satisfied for the present clusters, which are relatively dis- 
tant (Sect. [5]) and have been observed by 2MASS in a single 
epoch. 
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Figure 5. Same as Fig. |4] for the region R < 5' of LK 10. The 
equal-area comparison field extraction was taken from the re- 
gion 29.58' < R < 30'. A relatively populous giant clump and 
about 3.5 mag of the MS show up, especially in the decontami- 
nated CMDs, denoting advanced age. The 1 Gyr Solar-metallicity 
Padova isochrone is applied to the CMDs. Reddening vectors as 
in Fig. H 



4.1.1 Description of the algorithm 

We work with the st atistical algorithm introduced by 
iBonatto fc Bical (l2007al ) to deal with the field-star contami- 
nation in CMDs. The algorithm takes into account simulta- 
neously star-count density and colour/magnitude similarity 
between cluster and comparison field. It measures the rel- 
ative number densities of probable field and cluster stars 
in cubic CMD cells whose axes correspond to the J mag- 
nitude and the {J — H) and (J — Ka) colourfl The algo- 
rithm: (i) divides the full range of magnitude and colours 
covered by the CMD into a 3D grid, (ii) calculates the ex- 
pected number density of field stars in each cell based on 
the number of comparison field stars with similar magni- 
tude and colours as those in the cell, and (Hi) subtracts the 
expected number of field stars from each cell. By construc- 
tion, t he algorithm is sensitiv e to local field-star contami- 
nation (|Bonatto fc BicalbOOTal ). Typical cell dimensions are 
AJ = 1.0, and A(J - H) = A(J - Ks) = 0.25, which are 
large enough to allow sufficient star-count statistics in indi- 



^ These are the 2MASS colours that provide the maximum 
variance among CMD sequen ces for OCs of different ages (e.g. 



iBonatto. Bica fc Girardill2004) . 
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Figure 6. Same as Fig.|4]for the region _R < 4' of FSR 1521, with 
the equal-area comparison field extraction taken from 29.73' < 
R < 30'. The decontaminated CMDs suggest a giant clump and 
a red giant branch of an old distant OC. The age solution corre- 
sponds to the 2 Gyr isochrone. Reddening vectors as in Fig. |4] 



vidual cells and small enough to preserve the CMD morphol- 
ogy. The comparison fields are located within R = 20' — 40' 
(LKl), ii = 30'-60' (LKIO), and _R = 15' - 30' (FSR 1521 
and FSR 1555). The inner boundary of the comparison field 
lies beyond the probable tidal radius (Sect.[S|, which min- 
imises the probability of oversubtraction of cluster stars. We 
emphasise that the equal-area field extractions (middle pan- 
els of Figs. 13] -0 should be taken only for qualitative com- 
parisons. The decontamination is based on the large sur- 
rounding area as des cribed above. Fu rther details on the 
algorithm are given in iBonatto fc Bical (|2007a). 

As discussed in iBonatto fc Bical (|2007ai '). differential 
reddening between cluster and field stars may be critical for 
the decontamination algorithm. Important gradients would 
require large cell sizes or, in extreme cases, preclude appli- 
cation of the algorithm altogether. Basically, it would be 
required colour differences of \A{J — H)\ > cell size (0.25, 
in the present work) between cluster and comparison field 
for the differential reddening to affect the subtraction in a 
given cell. However, the cluster and comparison field CMDs 
(Figs. |4][7]) indicate that differential reddening is not impor- 
tant for the present sample. 




0.2 0.5 0.8 1.1 1.4 0.2 0.5 0.8 1.1 1.4 1.7 2.0 

(J-H) (J-K,) 



Figure 7. Same as Fig.|4]for the region R < 3' of FSR 1555. The 
equal-area comparison field extraction was taken from the region 
29.85' < i? < 30'. The decontaminated CMDs are best fitted 
with the 1.5 Gyr Solar-metallicity Padova isochrone. Reddening 
vectors as in Fig. l4l 



4.1-2 Decontamination statistics 

Additional statistical analysis is required because of the rela- 
tively high reddening values (Table [2} affecting the clusters. 
In Table [3] we present the full statistics of the decontami- 
nation, discriminated by magnitude bins. Statistically rele- 
vant parameters are: (i) Nia which, for a given magnitude 
bin, corresponds to the ratio of the decontaminated number 
of stars to the la Poisson fluctuation of the number of ob- 
served stars, (a) (TPS, which is related to the probability that 
the decontaminated stars result from the normal star count 
fluctuation in the comparison held and, (Hi) FSunii, which 
measures the star-count uniformity of the comparison field. 
Properties of A'"io-, aps, and FS uuis, measured in PCs and 
field f luctuations are discussed in lBica. Bonatto fc Camargol 
l|2008h . Table [3] also provides integrated values of the above 
parameters, which correspond to the full magnitude range 
spanned by the CMD of each OC. The spatial regions are 
those sampled by the CMDs shown in the top panels of 
Figs. HZl 

CMDs of star clusters shoul d have integrated Ni^ val- 
ues significantly higher than 1 (jBica. Bonatto fc Camargg 
I2OO8), a condition that is met for the present objects {Nic = 
6.1 — 8.9). As a further test of the statistical significance of 
the above results, we investigate star count properties of 
the field stars. First, the comparison field is divided into 8 
sectors around the cluster centre. Next, we compute the pa- 
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Table 3. Field-star decontamination statistics 
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Table Notes. For each magnitude bin (A J), we give the number of observed stars (A^obs) within the spatial region sampled in the CMDs 
shown in the top panels of Figs. l4l and [5l the respective number of probable member stars (A^cl) computed by the decontamination 
algorithm, the JVi^ parameter, the la Poisson fluctuation (o-ps) around the mean, with respect to the star counts measured in the 8 
sectors of the comparison field, and the field-star uniformity parameter. The statistical significance of A^cl is reflected in its ratio with the 
1(7 Poisson fluctuation of N^hs (A'^lo-) and with o-pg. The bottom line corresponds to the full magnitude range. The subtraction efficiency 
(fsub) is also given. 

rameter ups , which is the 1 a Poisson fluctuation around the 
mean of the star counts measured in the 8 sectors (corrected 
for the different areas of the sectors and cluster extraction). 
In a spatially uniform comparison field, crps is expected to 
be very small. Thus, OCs should have the probable num- 
ber of member stars (A'ci) higher than ~ Scrps, to minimise 
the probability that A*'ci arises from fluctuations of a non- 
uniform comparison field. This condition is fully satisfied by 
the present clusters, reaching the level A^ci ~ (5 — 10) upg. 
We also provide in Table|3]the parameter -F'5'unif . For a given 
magnitude bin we flrst compute the average number of stars 
over all sectors (A'^) and the corresponding la fluctuation 
cr(jv); thus, FSunif is deflned as FSunit = <^(n)/{N). Non 
uniformities such as heavy differential reddening should re- 
sult in high values of FSunii- On the other hand, FSunU is 
low (< 0.09), except for LKIO {FSnnii = 0.21), which re- 
flects a rather irregular fleld (Fig. [2]). Finally, we note that 
the number of decontaminated stars in most magnitude bins 
is larger than what could be expected from fleld-star fluctu- 
ations (A'ci > 3crps). 

Since we usually work with comparison fields larger 
than the possible-cluster extractions, the correction for the 
different spatial areas between field and cluster is expected 
to produce a fractional number of probable field stars (n^^") 



in some cells. Before the cell-by-cell subtraction, the frac- 
tional numbers are rounded off to the nearest integer, but 
limited to the number of observed stars in each cell n'^'^l = 
NI{nf") < nil", where NI represents rounding off to the 
nearest integer). The global effect is quantified by means of 
the difference between the expected number of fleld stars in 
each cell irffs^) and the actual number of subtracted stars 
{nTub)- Summed over all cells, this quantity provides an es- 
timate of the total subtraction efficiency of the process, 

fsub = 100 X ^ nltt I n)f (%) . 

cell cell 

Ideally, the best results would be obtained for an efficiency 
faub ~ 100%. With the assumed grid settings for the decon- 
tamination of the present clusters, the subtraction efficien- 
cies turned out to be higher than 90%. 

4-1.3 Decontaminated CMDs and surface densities 

As an indicator of the algorithm efficiency we can take the 
decontaminated stellar surface-density distributions (bot- 
tom panels of Figs. [5] and [3}. The central excesses have 
been significantly enhanced with respect to the raw photom- 
etry (top panels), while the residual surface-density around 
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the centre has been reduced to a minimum level. By de- 
sign, the decontamination depends essentially on the colour- 
magnitude distribution of stars located in different spatial 
regions. The fact that the decontaminated surface-density 
presents a conspicuous excess only at the assumed cluster 
position implies significant differences among this region and 
the comparison field, both in terms of colour-magnitude and 
number of stars within the corresponding colour-magnitude 
bins. This meets the expectations of star clusters, which 
can be characterised by a single-stellar population, projected 
against a Galactic stellar field. 

The decontaminated CMDs are shown in the bottom 
panels of Figs. |4]-[7l As expected, essentially all of the disc 
contamination is removed, leaving stellar sequences typical 
of reddened old OCs, with well-developed red clumps and 
different extents of the MS. 

As a caveat, we cannot exclude the possibility of differ- 
ential reddening to account for part of the observed spread 
in the CMD distribution of stars. To examine this issue we 
show, in the bottom panels of Figs. |4][7l reddening vectors 
computed with the 2MASS ratios (Sect. 2]) for a standard 
visual absorption Ay = 2. Given the absorptions derived for 
the clusters (Table [2]), the standard value, which represents 
from ~ 25% to « 50% of the total Ay, can be taken as 
an upper limit to the differential reddening. Together with 
the decontaminated CMDs, this experiment suggests that 
differential reddening in all cases is not large, because the 
giant clumps are rather tight, while the MS spread appears 
to be dominated by photometric errors. The typical depen- 
dence of the 2MASS photometric errors on magnitude, for 
objects projecte d along different directions, is discussed in 
iBonatto fc Bical (|2007b). 

We conclude that the qualitative and quantitative ex- 
pectations of the decontamination algorithm have been sat- 
isfied by the output. In all cases, the decontaminated pho- 
tometry presents a conspicuous excess, with respect to the 
surroundings, in the surface-density distribution (Figs. [2] 
and [3]). In addition, field-decontaminated CMDs extracted 
from the spatial regions where the excesses occur (Figs.|3]lZl)j 
present statistically significant (Table [Sj cluster CMDs. 



to preserve the several peculiarities present in the TP-AGB 
tracks, namely, the cool tails of C-type stars due to the use 
of proper molecular opacities as convective dredge-up oc- 
curs along the TP-AGB, the bell-shaped sequences in the 
HR diagram for stars with hot-bottom burning, the changes 
of pulsation mode between fundamental and first overtone, 
the sudden changes of mean mass-loss rates as the surface 
chemistry changes from M- to C-type, etc. Because it is im- 
portant to derive, or at least set constraints on astrophysical 
parameters of old star clusters, we adopt as working strategy 
the search for solutions wi thin a range of ages and metallic- 
ities. As discussed by, e.g. iFriell |l995), OC metallicities in 
general range from Solar {[Fe/ H] = 0, or Z = 0.019) to sub- 
Solar {[Fe/H] = -0.5, Z = 0.006, or ^ 1/3 S olar) value s. 
According to the metallicity gradient (Fig. 7 in lFriel|[l995l ). 
the present OCs fall around the locus occupied by OCs with 
metallicity [Fe/H] ~ —0.1. Thus we base the following anal- 
ysis on this most probable value. To com pute Galactocen tric 
distances, we adopt Rq = 7.2±0.3kpc l|Bica et al.ll2006l ) as 
the Sun's distance to the Galactic centre. This value was 
derived by means of the GC spatial distributior(f|. 

Historically, different approaches have been used to ex- 
tract astrophysical parameters from CMDs by means of 
isochrone fits. The simplest ones are based on a direct com- 
parison of a set of isochrones with the CMD morphology, 
while the more sophisticated include photometric uncer- 
tainties, binarism, and variations on metallicity. Most o f 
these methods are summarised in iNavlor fc Jeffries! (|2006l ). 
in which a maximum-likelihood CMD fit method is de- 
scribed. 

For simplicity, in the present cases fits are made by eye, 
with the tight giant clumps as the strongest constraint. We 
also require that, because of the potential presence of bi- 
naries, the adopted isochrone should be shifted somewhat 
to the left of the MS fiducial line, i.e. a median line that 
takes into account the MS spread, including the photometric 
uncertainties as well (e.g. iBonatto. Bica fc Santos Jr. I l2005l . 
and references therein). In the following Section we discuss 
each cluster individually. 



5 CLUSTER AGE, REDDENING, AND 
DISTANCE 

The field-decontaminated CMD morphologies derived in 
Sect. U can be used to compute the cluster funda- 
mental parameters. We work with Padova isochrones 
(|Girardi et al.l \200^ ) computed with the 2MASS 
J, H, and Kg filter^. The updated isochrones 
are very similar to the Johnson-Kron-Cousins ones (e.g. 
iBessel fc Brett 1988h. with differences o f at most 0.01 in 
(J - H) (|Bonatto. Bica fc GirardU bOO^ ). Distinctive fea- 
tures of the updated isochrone set are centred mostly on 
the greatly-improved treatment of the thermally-pulsing 
asymptotic g i ant b ranch (TP-AGB) phase. According to 
iMarigo et al.l l|2008l ). the updated isochrones are intended 



5.1 FSR1521 

The decontaminated CMD morphology of FSR1521 (Fig.[6} 
shows a concentration of stars at J « 15, which indicates a 
relatively populous MSTO of an old cluster detected near 
the 2MASS photometric limit. Besides, a distinctive red 
clump shows up at J « 12.5 and (J — H) « 0.9. Taken 
together, these stellar sequences characterise a distant old 
OC. 

Allowing as well for the photometric uncertainties, ac- 
ceptable fits to the decontaminated CMD morphology are 
obtained with the 2 Gyr isochrone, with an uncertainty of 
±0.5 Gyr. This can be seen in Fig. |8] in which we com- 
pare solutions with the ages 1.5, 2, and 2.5 Gyr, both for 
Solar and sub-Solar {[Fe/H] — —0.5) metallicities. Ages 
younger than about 1.5 Gyr would produce a significantly 



\http://stev.oapd. inaf.it/cgi-bin/cmd\ - Bolometric and colour 
corrections were computed for a set of isochrones using the 
2MASS filter responses. The isochrones were subsequently pro- 
vided in the Vega Mag system. 



^ Other recent studies g ave similar results, e.g. Rp, = 7.2ili0.9kpc 
IIEis enhauer et al ] |200ai Rq = 7.6 2 ill 0.32 kpc j Eisenhauer et al.l 
2005^ and Rq = 7.52 ± O.lOkpc llNishivama~t al.. .20061) . with 
different approaches. 
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Figure 8. Possible solutions for the age and metallicity of 
FSR1521. Padova isochrones with the ages 1.5 (left panels), 2 
(middle), and 2.5 Gyr (right) are used. Except for the Solar metal- 
licity 1.5 Gyr isochrone, the remaining ones provide acceptable 
fits, either with Solar- metallicity (top panels) or the 1/3 Solar- 
metallicity (bottom) solutions. 



poorer description (much fainter than observed) of the red 
clump, while those older than about 2.5 Gyr, would be pro- 
hibitively shifted to the right of the MSTO. In any case, 
both metallicity ranges provide similar fits to the decon- 
taminated morphology. Since the remaining OCs have ages 
(Table [2]) and CMDs with photometric uncertainties sim- 
ilar to those of FSR1521, the above test shows that any 
metallicity within 1/3 < Z/Zq < 1 would produce accept- 
able solutions. Thus, for simplicity and because the present 
OCs are relatively near the Solar circle, we adopt the Solar- 
metallicity isochrones as more probable solutions. 

We caution that, because of the 2MASS photometric 
uncertainties for the lower sequences, a more sophisticated 
approach for isochrone fitting might produce an overinter- 
pretation. For this reason, we apply the above method for 
all objects. 

With the adopted solution, the fundamental parameters 
of FSR1521 are the near-infrared reddening E{J — H) = 
0.34 ± 0.05, which corresponds to E{B -V) = 1.1 ± 0.2, or 
Ay = 3.4 ± 0.5, the observed and absolute distance moduli 
(m - M)j = 14.2 ± 0.3 and (m - M)o = 13.26 ± 0.33, re- 
spectively, and the distance from the Sun = 4.5±0.7kpc. 
Thus, for _R0 = 7.2 kpc, the Galactocentric distance of 
FSR1521 is Rgc = 7.8 ± 0.4 kpc, which puts it ^ 0.6 kpc 



outside the Solar circle. This solution is shown in Fig. [6] 
(bottom panels). 

We call attention that the best fit takes into account 
strong constraints provided by the tight distribution of the 
clump giants, which preserves accuracy especially for the 
reddening and distance. Thus, the resulting uncertainties are 
adopted as the most probable errors, which are subsequently 
propagated to the different sets of parameters (Tables [2] and 
|4|. However, given the uncertainties in the faint-photometric 
end, age is less constrained. 



5.2 LKl 

The best-fit to the CMD of LK 1 was obtained with the 
1 Gyr isochrone, E{J ~ H) = 0.90 ± 0.03 and (m - M)j = 
15.5 ± 0.2. Taking into account fit uncertainties we derive 
the age 1.0±0.2Gyr, E{B - V) = 2.9±0.1, Ay = 8.9±0.3, 
(m - M)o = 13.02 ± 0.22, do = 4.0 ± 0.4 kpc, and Rgc = 

7.3 ± 0.2 kpc. LKl is « 0.1 kpc outside the Solar circle. 
This solution is shown in the bottom panels of Fig. |4l It 
represents well the MSTO, red clump and giant branch, in 
both colours. The few stars bluer than the MSTO might be 
blue stragglers or un-subtracted field stars. Our derived dis- 
tance modulus is about 2 mags f ainter than that adopted by 
iLe Duigou fc Knodlsede? (200?). 

With 1-2 kpc as the distance to the Cygnus Association 



Le Duigou fc Knodlsedei1l2002l : lBica. Bonatto fc Dutral 



l2003h . LKl, at « 4kpc, is a distant background old open 



cluster in the line of sight of that association. 



5.3 LKIO 

The decontaminated CMD morphology of LK 10 (Fig. [5} 
is somewhat more constrained than that of LKl. Besides 
a tight red clump, it features about a 3 mag MS extent. 
Fundamental parameters of LK 10 are an age 1.0 ± 0.1 Gyr, 
E{J - H)= 0.86 ± 0.02, which corresponds to E{B - V) = 
2.8 ± 0.1 and Av = 8.5 ± 0.3, (m - M)j = 13.1 ± 0.1, 
(m - M)o = 10.7 ± 0.11, dQ = 1.4 ± 0.1 kpc, and Rgc = 
7.1 ± 0.1 kpc. LKIO lies « 0.1 kpc inside the Solar circle. 
Our derived distance modulus a grees with that adopted by 
iLe Duigou fc Knodlsedeil (|2002l ). 

With ~ 1.4 kpc, LKIO might be spatially coinci- 
dent with the Cygnus Association. However, they must not 
share a common origin, since LK 10 is much older than 
the estimated age of the Cygnus Association , 1-4 Myr (e.g. 
iMassev. Johnson fc DeGioia-EastwoodI Il995l ). Despite the 
different distances, LK 1 and LK 10 are affected essentially 
by the same reddening value. 



5.4 FSR1555 

With the MSTO and red clump detected in the decon- 
taminated CMD of FSR1555 (Fig. 0, we derive an age 

1.5 ± 0.5 Gyr, E{J - H) = 0.37 ± 0.02, which corresponds 
to E{B -V) ^ 1.2 ± 0.1 and Av = 3.7 ± 0.3, (m - M) j = 
14.1 ± 0.1, (m - M)o = 13.08 ± 0.11, dg = 4.1 ± 0.2kpc, 
and Rgc = 7.1 ± 0.1 kpc. FSR 1555 lies « 0.1 kpc inside the 
Solar circle. 
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6 CLUSTER STRUCTURE 

Structural parameters are derived by means of the projected 
radial density profiles (RDP) built with the stellar num- 
ber density around the cluster centre. Usually, star clus- 
ters have RDPs that follow a well-defined analytical pro- 
file. Among these ar e the empiri cal, single mass, modified 
isotherm al sphere of iKind (|f 9661 ). the modified isothermal 
sphere of lWilsonI l| 19751 ). which assumes a pre-defined stel- 
lar distribution function (and produces more extended en- 
ve lopes than King 1966) , and the power law with a core 
of[El son. Fall fc Freemanl l| 19871 ). These functions are char- 
acterised by different parameters that are related to clus- 
ter structure. However, considering the error bars of the 
present RDPs (Fig. [9]), we adopt the analytical function 
a(i?) = abg -|- cro/(l + [R/RcY), where Gbg is the residual 
background density, ao is the central density of stars, and 
Rc is the co re radius. T his function is similar to that in- 
troduced by iKind (| 19621 ) to describe the surface brightness 
pr ofiles in the central parts of globular clusters. As discussed 
in iBonatto fc Bical (|2008al ). RDPs built with depth- limited 
photometry produce structural radii comparable to the in- 
trinsic (i.e. derived with deep photometry) ones. 

To minimise noise in the RDPs, we first apply a colour- 
magnitude filter to the photometry, which excludes stars 
with colours unlike those of the cluster sequence. Colour- 
magnitude filters are wide enough to include cluster MS and 
evolved star colour distributions, as well as the la photo- 
metric uncertaintie^. The colour-magnitude filters for the 
present OCs are shown in the bottom-left panels of Figs. [4]- 
[71 However, residual field stars with colours similar to those 
of the cluster are expected to remain inside the colour- 
magnitude filter. They affect the intrinsic stellar RDP in a 
way that depends on the relative densities of field and clus- 
ter stars. The contribution of the residual contamination 
to the observed RDP is statistically subtracted by means 
of the field. As a result, the use of colour-magnitude filters 
enhances the contrast of the RDP with res pect to the back- 
groun d, especially in crowded fields (e.g. iBonatto &: Bical 
l2007al ). 

Oversampling near the centre and undersampling at 
large radii are avoided by using rings of increasing width 
with distance from the cluster centre. A typical set of 
ring widths is A_R = 0.5, 1, 2.5, and 5', respectively for 
0' «; < 1', 1' s; i? < 4', 4' s; i? < lO', and R > 10'. The 
number and width of the rings can be set to produce RDPs 
with adequate spatial resolution and small la Poisson errors. 
The residual background level of each RDP corresponds to 
the average number of colour-magnitude filtered stars mea- 
sured in the field. The R coordinate (and uncertainty) of 
each ring corresponds to the average position and standard 
deviation of the stars inside the ring. 

The colour-magnitude filtered RDPs of the present clus- 
ters are shown in Fig. |9l where we also show the pro- 
files produced with the observed (raw) photometry. As ex- 



^ Colour-magnitude filter widths should also account for forma- 
tion or dynamical evolution-related effects, such as enhanced frac- 
tions of binaries (and other multiple systems) towards the central 
parts of c lusters, since such systems tend to wide n the MS (e.g . 



pected, minimisation of the number of non-cluster stars by 
the colour-magnitude filter resulted in RDPs with higher 
contrast with respect to the background. Fits of the King- 
like profile were performed with a non-linear least-squares 
fit routine that uses errors as weights. To minimise degrees 
of freedom, uo and Rc were derived from the RDP fit, while 
abg is measured in the field. The best-fit solutions are shown 
in Fig. [9] and the fit parameters are given in Table|4l For ab- 
solute comparison with other clusters. Table 2] also presents 
parameters in absolute units, based on the cluster distances 
(Sect. [5]). Because of the 2MASS photometric limit, which 
for the present clusters corresponds to a cutoff for stars 
brighter than J ~ 16.5, cto should be taken as a lower limit. 

Within uncertainties, the adopted King-like function 
describes well the colour-magnitude filtered RDPs along 
the full radius range, especially for LKl, FSR1521, and 
FSR 1555. The exception is LK 10, which shows a marked ex- 
cess in the central region. This central cusp in LK 10 suggests 
a post-core collapse phase in this r-^ 1 Gyr OC, as detected in 
part of the globular clusters fe.g. lTrager. King fc Diorgovskil 
I1995I ). Such a central cusp has been observed in the RDP of 
other OCs, such as NGC3960 tBonatto fc B ica 2006). 

We also estimate the cluster radius (-Rrdp) by visu- 
ally comparing the RDP level (and fluctuations) with the 
background. It corresponds to the distance from the clus- 
ter centre where RD P and ba ckground are statistically in- 
distinguishable (e.g. IBonatto fc Bica 2005,. and references 
therein). Thus, most of the cluster stars are contained within 
^?RDP, which should not be mistaken for the tidal radius. 
Tidal radii are derived from, e.g. the 3-parameter King- 
profile fit to RDPs (see below), which requires large sur- 
rounding fields and adequate Poisson errors. For instance, 
in populous and relatively high Galactic latitude OCs such 
as M 67, NGC 188, and NGC 2477, cluster radii are a fac- 
tor ~ 0.5 — 0.7 of the respective tidal radii (IBonatto fc Bical 
iooi). The factor is somewhat lower for low-Galactic lati- 
tude objects, because of the enhanced background. The clus- 
ter radii of the present objects are given in cols. 5 (angular 
scale) and 12 (absolute scale) of Table [l] 

Table [4] (col. 7) provides the density contrast parame- 
ter 5c = 1 -\- ao/abg, which, for the present clusters is rela- 
tively high (3.3 < 5c < 8.5). Since 5c is measured in colour- 
magnitude-filtered (lower noise) RDPs, it is usually higher 
than the visual contrast produced by images (e.g. Fig. [TJ. 

Alternatively, we tried t o fit the R DPs with the 3- 
parameter function (based on iKindl 19621 ) 



a{R) = (Jo 



^1 + [R/Rcy yjl + {RtlRcY 



Hurley fc T out 1998; Kcrbcr ct al. 2002. ; .Bonatto fc Bicall2007al : 



Bonatto. Bica fc Santos Jr.ll2005h . 



which includes the tidal radius (-Rt)- However, convergence 
occurred only for LKl, with Rt = 18.3' ± 5.5'. With 
such radii, the concentration parameter of LK 1 is c = 
log(_Rt/-Rc) ~ 1-5 which, compared to the Galactic (non-core 
colla pse) globular clusters, p uts it around the median value 
(e.g. IBonatto fc Bicall2008bl ) . Qualitatively, the 3-parameter 
fit to the RDP of LK 1 is indistinguishable to that shown in 
Fig. |9l within the uncertainties. Thus, we simply show in 
panel (c) the position of the tidal radius. 

Interestingly, the cluster radii of LK 1 and LK 10 given 
bv lLe Duigou fc Knodlsedeil (2002) arc about half those de- 
rived in this work. This probably occurs because we work 
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Table 4. Derived cluster structural parameters 



Cluster 






Rc 


Rrdp 


Rt 


Sc 


1' 




O"0 


Rc 


Rrdp 


Rt 






{*'-') 


(') 


(') 


(') 




iPc) 


(*pC-2) 


(*pc~2) 


(pc) 


(pc) 


(pc) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


LKl 


5.15 ±0.03 


38.5 ±8.3 


0.57 ±0.10 


6.0 ±1.0 


18 ± 5 


8.5 ±1.6 


1.163 


3.8 ±0.1 


28.5 ±6.1 


0.66 ± 0.12 


7.0 ±1.2 


21 ±6 


LKIO 


2.67 ±0.02 


9.5 ± 3.3 


1.94 ±0.57 


16.0 ±2.0 




4.6 ±1.2 


0.444 


16.3 ±0.1 


58.1 ± 17.1 


0.78 ± 0.23 


6.5 ±0.8 




FSR1521 


6.66 ±0.06 


15.4 ±7.1 


0.64 ±0.24 


5.2 ±1.0 




3.3 ±1.0 


1.302 


3.9 ±0.1 


9.1 ±4.2 


0.83 ± 0.31 


6.9 ±1.3 




FSR1555 


7.62 ± 0.03 


37.8 ±2.6 


0.53 ±0.05 


4.5 ±0.5 




6.0 ±0.3 


1.197 


5.3 ±0.1 


26.4 ± 1.8 


0.63 ± 0.06 


5.4 ±0.6 





Table Notes. Core (Rc), cluster (-Rrdp), a^nd tidal (Rt) radii are given in angular and absolute units. Col. 7: cluster/background density 
contrast parameter {Sc = 1 + fo/fbg), measured in the colour-magnitude filtered RDPs. Col. 8: arcmin to parsec scale. 
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Figure 9. Stellar RDPs built with colour-magnitude filtered 
(filled circles) and raw (empty) photometry. Solid line: best-fit 
King-like profile. Horizontal shaded polygon: background stellar 
level. Shaded regions: la King fit uncertainty. The core {Rc) and 
cluster (i?RDp) radii are indicated in all cases, while for LK 1, the 
tidal radius (-Rt) is also shown. Note the central density excess in 
LK 10 profile. 



with colour-magnitude-filtered RDPs, which enhances the 
cluster/background contrast and probes larger cluster ex- 
tensions. With respect to FSR 1521 and FSR 1555, our val- 
ues for the core radius are abou t 1/3 of those given by 
iFroebrichT Scholz fc Raftervl l|2007h . while our -Rrdp is con- 
sistent with their Rt , given the above relation between both 
radii. 

Compared to the distribution of core radius der i ved fo r 
a sample of relatively nearby OCs bv lPiskunov et all ((20o3), 
the present OCs, especially FSR 1521, occupy the small- -Rc 



tail. Besides, for a relation between tidal and cluster radius 
as ~ 2 X Rnup, LK 10, FSR 1521, and FSR 155 are located 
around the median value, while LK 1 occupies the large- -Rt 
tail. 



7 MASS ESTIMATES 

Since we detect ~ 3 mags of the MS of LK 10, we can 
build its mass function (0(m) = and compute the 

mass stored in stars. We work with colour-magnitude fil- 
tered photometry, the 3 2MASS bands separately, and 
the mass-luminosity relation obtained from the corre- 
sponding Padova isochrone and distance from the Sun 
(Sect. [5]). Fu rther details on MF construction are given in 
iBica. Bonatto fc Blumbergj l|2006l ). The effective MS stellar 
mass range is 1.15 ^ m{MQ) ^ 1.80. 

The MF of LK 10 is well represented by the function 



m) oc m 



-(i+x) 



with the slope x = 2.4 ± 0.4. Within 
un certainties, this slope agrees wit h x = 1-98 ± 0.25 derived 
bv lLe Duigou fc Knodlsedeij (|2002|). Both values are steeper 
than the y = 1.35 of ISalpetej (|l955l ) initial mass function 
(IMF). 

The number of observed MS and evolved stars in LK 10 
(for -R ^ -Rrdp) is derived by counting the stars (in the 
background-subtracted colour- magnitude filtered photome- 
try) that are present in the mag ranges, 13.4 < J < 16.4 
for the MS and J < 13.4 for the evolved stars. There are 
nMs ~ 959 ± 20 and Ucvoi = 25 ± 8, MS and evolved stars, 
respectively; the corresponding mass values are niMS = 
1311 ± 26 Mq and nicvoi = 46 ± I6M0. The evolved star 
mass corresponds to ricvoi multiplied by the stellar mass at 
the TO, rriTO = 1.8 Mq. Thus, the observed stellar mass of 
LK 10 is iriobs ~ 13 60 -Mp), which agrees with the eq uivalent 
value estimated bv lLe Duigou fc Knodlsederl (|2002t ). 

Finally, we estimate the total stellar mass by extrapo- 
lating the observed MF down to the H-burning mass limit 
(0.08 Mq). We follow the universal IMF of Kroupa (200l]), 
which assumes increasing fiattening towards low-mass stars. 
This IMF is described by the slopes x = 0.3 ± 0.5 for the 
range 0.08 < m(MQ) 0.5 and x = 1-3 ± 0.3 for 0.5 < 
m(Af0) ^ 1.0. We obtain ruc^tr = 4420 ± 2000 Mq. Thus, 
the stellar mass of LK 10 can be put in the range 1360 — 
4420 Mq. Again, within uncertai nties, the upper value is 
consistent with that estimated bv lLe Duigou fc Knodlsederl 
(2002). 

Because of the limited MS range, only estimates of the 
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observed cluster mass are made for the remaining objects, 
by means of the age solutions given in Sect. (5] In all cases we 
consider the region within R -Rrdp (Table |4]). For LK 1 
we derive rriMS = 108 ± 9 Mq and rrievoi = 270 ± 26 M©, 
which leads to the total observed mass mobs ~ 380 Mq . This 
value corresponds to about 1/ 3 of the estimate given by 
iLe Duigou fc Knodlsedeil l|2002l ). The values for FSR1555 
are tums = 123 ± 9 Mq, nievoi = 139 ± 16 Mq, and niobs ~ 
260 Mq. Since FSR 1555 and LK 1 are approximately at the 
same distance from the Sun, the latter appears to be some- 
what more massive than the former. 

With the more probable solution (2Gyr) for FSR 1521, 
the total number of stars present in ii < 5.2' CMD is 
ntot ~ 161, of which « 18 are in the red clump and ~ 111 
in the MSTO. Thus, for a MSTO mass of « 1.63 M©, the 
observed mass of FSR 1521 is mobs ~ 262 Mq, of which 
« 81 Mq correspond to the evolved stars, and ~ 29 Mq are 
stored in the red clump. Interestingly, while the observed 
mass of FSR 1521 is similar to that of FSR 1555, the evolved 
mass is somewhat lower. Since both OCs are at comparable 
distances, this difference is consistent with the relative ages. 

The fact that LK 1, FSR 1521, and FSR 155 present sim- 
ilar observed masses as LK 10, suggest that they all might 
be as massive as the latter OC. 

Given the accuracy of the isochrone best-fit for old OCs 
provided by near-infrared decontaminated CMDs, e.g. Fig. [8] 
and related discussions, the observed masses are expected to 
be representative. Of course, deeper near-infrared photom- 
etry coupled to similar methods as the present one, and/or 
including proper motion filtering, would produce more con- 
strained results. 



8 DISCUSSION 

With the analyses of the preceding sections we have gath- 
ered important clues to establish that the objects dealt with 
in this paper are Gyr-class OCs, or older. We also derived 
representative, i.e. constrained by means of the isochrone fit 
(Sect. [5|, fundamental and structural parameters, most of 
which have been derived for the first time. We use these pa- 
rameters to put the clusters into perspective, by comparing 
some of their properties with those of a set of well-studied 
OCs. 

As reference we take the nearby OCs with ages in the 
range 70 - 7000Myr and rn a sses w ithin 400 - 5 3OOM0 
studied by iBonatto Sz Bica ll2005h, together with the 
young OCs N GC6611 llBonatto. Santos Jr. fc Bical |2006| ) 
and NGC4755 l|Bonatto et al.ll2006bl 'l. The reference clus- 
ters are distinguished according to tot a l mas s (higher or 
lower than 1 000 Mq). iBonatto fc Bical (|2005h discuss pa- 
rameter correlations in the reference sample. 

As shown in panel (a) of Fig. 1101 the core and clus- 
ter radii of the reference OCs are related by 7?rdp = 
(8.9 ± 0.3) X i?Le*° '', which suggests a similar scaling 
for both kinds of radii, at least for the radii ranges 0.2 < 
■Rc(pc) < 1.5 and 2 < i?RDp(pc) < 15. LKl, LKIO, 
FSR 1521, and FSR 1555 fit tightly in the relation. They 
also appear to follow the trend of increasing cluster radii 
with Galactocentric distance (panel b). T his kind ofde- 
pendence was previously suggested by, e.g. iLvngal l|l982l ). 
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Part of this relation may be primordial, in the sense that 



Figure 10. Relations involving OC structural and fundamen- 
tal parameters. Circles: reference OCs. Dotted circles: massive 
(> 1 000 Mq) OCs. We caution that LKIO may be a post-core 
collapse open cluster. In this figure we adopted the core radius 
derived from the King-like fit. 



the higher molecular gas density in central Galactic regions 
m ay have produced clusters with sm a ll rad ii, as suggested 
bv Ivan den Bergh. Morbev fc Pazderl (|l99lh to explain the 
increase of globular cluster radii with Galactocentric dis- 
tance. After formation, mass loss associated with stellar and 
dynamical evolution (such as mass segregation and evapo- 
ration), together with tidal interactions with the Galactic 
potential and giant molecular clouds, also contribute to the 
depletion of star clusters, especially the low-mass and cen- 
trally located ones (Sect.[l|. A similar dependence on Galac- 
tocentric distance for Rc is implied by the data shown in 
panel (a). 

In panels (c) and (d) of Fig.[Tn]we compare the presently 
derived cluster and core radii with those of the refer- 
ence sample in terms of age. LKl, LKIO, FSR 1521, and 
FSR 1555 have core and cluster radii similar to those mea- 
sured in the reference OCs of equivalent age. 

We show in Fig. [TT] the spatial distribution of LKl, 
LK 10, FSR 1521, and FSR 1555, as they lay in the Galactic 
pl ane. The spiral arm st ruct ure of the Milky Way is ba sed 
on lMomanv et al.l (120061 ) and lDrimmel fc Spergell (l200lh. as 
deriv ed from HII regions, and molecular clouds fe.g. iRusseill 
[200t ). The Galactic bar is shown with an orientation of 14° 
and 6kpc in total length (lFreudenreichlll998l : IValleell2005l ). 
The present OCs are compared to the spatial distribution 
of the OCs with known age given in the WEBDA database. 
For comparison purposes we consider two age groups, clus- 
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Figure 11. Spatial distribution of the present star clusters com- 
pared to the WEBDA OCs with ages younger (gray circles) 
and older than 1 Gyr (black dots) . Clusters are overplotted on 
a schematic projection of the Galaxy, as seen from the North 
pole, with 7.2 kpc as the Sun's distance to the Galactic centre. 
Main structures are identified. 



ters younger and older than 1 Gyr (the old OCs listed in 
Table [1] are merged into this group). As expected, old OCs 
are found preferentially outside the Solar circle, while the 
inner Galaxy contains few OCs so far detected. Besides, be- 
cause of the presence of bright stars, young OCs can be 
detected farther than the old ones, especially towar ds the 
central region. As discussed in lBonatto et al.] l|2006al ). cen- 
tral directions farther than « 2 kpc begin to be critically 
affected by incompleteness (due to crowding and high back- 
round levels) and enhanced disruption rates, while the drop 
in the number of OCs towards the opposite direction ap- 
pears to be mostly real. Also, all directions show depletion 
in the number of detected OCs farther than ~ 2 kpc. Consid- 
ering the overall distribution, the pairs LK 1 and LK 10, and 
FSR 1521 and FSR 1555 probe basically unexplored regions, 
in opposite directions. FSR 1555, and especially FSR 1521, 
are located close to the Carina-Sagittarius Arm. LK 1 and 
LK 10 are found between the Carina-Sagittarius and Orion- 
Cygnus arms. 



9 SUMMARY AND CONCLUSIONS 

Accurate fundamental and structural parameters of old star 
clusters are important for reasons that range from the com- 
pleteness of the open cluster parameter space to the deter- 
mination of the efficiency of cluster dissolution mechanisms. 
Taken together, the results of the present paper, and those 
of previous works by our group (Table [T|, add to a total 
of 18 old (age > 0.8 Gyr) open clusters with reliabl4i3 pa^- 



In the context discussed in Sect. [5] with uncertainties propa- 
gated from the best-fitting 2MASS isochrones. 



rameters derived, for the first time for most of them. New 
findings and the age-determination within a reasonable con- 
fidence level are important as well to improve the statistical 
coverage of the open cluster databases, especially for the def- 
inition of the old-age tail of the open cluster age-distribution 
function. 

In this paper we use colour-magnitude diagrams and 
radial density profiles to derive fundamental and structural 
parameters of the infrared open clusters LK 1 and LK 10, 
as well as of the star cluster candidates FSR 1521, and 
FSR 1555. Our approach is essentially based on field-star 
decontaminated 2MASS photometry, which enhances clus- 
ter CMD evolutionary sequences, and produces more con- 
strained parameters. 

We present consistent evidence, in the form of CMD 
morphology, statistical tests, structural parameters, and 
comparison with nearby OCs, that the objects are Gyr-class 
open clusters. With absorptions in the range 3.4 ^ Ay ^ 
8.9, the objects are well suited for the 2MASS photometry. 
They are located at ~ 4 kpc from the Sun, except for 
the relatively nearby LK 10, which is at ~ 1.4 kpc. The 4 
clusters are within ~ 0.6 kpc from the Solar circle. The old 
ages of LK 1 and LK 10 (and the distance of the former) pre- 
clude a common origin with the Cygnus Association. With 
a mass function slope (x = 2.4 ± 0.4) somewhat steeper 
than Salpeter's {% ~ 1.35), LK 10 is relatively massive, with 
a total mass within 1360 ^ m(MQ) ^ 4400. The remain- 
ing clusters present CMD morphology and observed mass 



(within 260 < m{M(7 



380) similar to those of LK 10, 



which suggests that they may be old clusters as massive as 
LK 10. This is consistent especially with the relatively well- 
populated red clumps. 

Structurally, LK 1, FSR 1521, and FSR 1555 are charac- 
terised by core and cluster radii similar to those of a sample 
of nearby open clusters of comparable age. This is consistent 
with the fact that, because they are relatively close to the 
Solar circle, they suffer similar tidal stresses. LK 10, on the 
other hand, presents evidence of post-core collapse, a struc- 
ture already observed in other old OCs, such as NGC 3960 
(Bonatto & BLca 2006,). 

Finally, we note that systematic searches in catalogues 
of star cluster candidates, coupled with an efficient field- 
star decontamination algorithm, can provide important ad- 
ditions the the known population of open clusters, old ones 
in particular. 
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